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Photons carrying non-zero orbital angular momentum (twisted photons) are well- 
Qh| known in optics. Recently, it was suggested to use Compton backscattering to boost 

optical twisted photons to high energies. Twisted electrons in the intermediate energy 
range have also been produced recently. Thus, collisions involving energetic twisted 
\ particles seem to be feasible and represent a new tool in high-energy physics. Here we 

J> ' discuss some generic features of scattering processes involving twisted particles in the 

initial and/or final state. In order to avoid additional complications arising from non- 
ly-j , trivial polarization states, we focus here on scalar fields only. We show that processes 

involving twisted particles allow one to perform a Fourier analysis of the plane wave 
cross section with respect to the azimuthal angles of the initial particles. In addition, 
using twisted states one can probe the autocorrelation function of the amplitude, which 
is inaccessible in the plane wave collisions. Finally, we discuss prospects for experimental 
study of these effects. 

x. 

g3 ■ 1 Introduction 

In perturbative quantum field theory we assume that interaction among the fields can be 
treated as a perturbation of the free field theory. This perturbation leads to scattering be- 
tween asymptotically free multiparticle states, which are usually constructed from the plane 
wave one-particle states. This choice greatly simplifies the calculations and represents a very 
accurate approximation to the real experimental situation in virtually all circumstances. How- 
ever, one can, in principle, choose any complete basis for the one-particle states other than the 
plane wave basis, provided that it is still made up of solutions of the free field equations. Such 
states can carry new quantum numbers absent in the plane wave choice and, if experimentally 
realized, they can offer new opportunities in high-energy physics. 

Thanks to the progress in optics made in the last two decades, it is now possible to create 
laser beams carrying non-zero orbital angular momentum (OAM) [lj, for a recent review see 
j2]. The lightfield in such beams is described via non-plane wave solutions of the Maxwell 
equations. Each photon in this lightfield, which we call a twisted photon, carries a non-zero 
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OAM quantized in units of h. Several sets of solutions have been investigated, such as Bessel 
beams or Gauss-Laguerre beams, but in all cases the spatial distribution of the lightfield 
is necessarily non-homogeneous in the sense that the equal phase fronts are not planes but 
helices. Such states form a complete basis which can be used to describe the initial and final 
asymptotically free states. Moreover, it is the basis of choice for experimental situations when 
the initial states are prepared in a state of (more or less) definite OAM. 

Twisted photons have been produced in various wavelength domains, from radiowave [3] 
to optical, with prospects to create a brilliant X-ray beam of twisted light in the keV range [4J . 
Very recently it was suggested to use the Compton backscattering of twisted optical photons 
off an ultra-relativistic electron beam to create a beam of high-energy photons with non-zero 
OAM [HIE]. The technology of Compton backscattering is well established [7j, and the high- 
energy electron beams and the OAM optical laser beams are already available. In addition, in 
the last months several groups have reported successful creation of twisted electrons, first using 
phase plates [8] and then with computer-generated holograms [i9j. Twisted electrons carried 
the energy as high as 300 keV and the orbital quantum number up to ~ 100. With all these 
achievements, creating high-energy particles in a controlled orbital angular momentum state 
and colliding them seems now feasible. It is therefore very timely to ask what new insights 
into the properties of particles and their interactions one can gain with this new degree of 
freedom. 

In this paper we begin this exploration by studying several generic scattering processes 
involving twisted particles in the initial and/or final states. Namely, we consider three specific 
cases: 

• single-twisted scattering: collision of a twisted state with a plane wave, 

• double-twisted scattering: collision of two twisted states, 

• two-particle decay of an unstable twisted particle. 

In the first two cases we assume that the final system X is described by plane waves, while 
in the last case we consider three choices for the final two-particle state: when both particles 
are plane waves, when one is twisted, and when both are twisted. The single- and double- 
twisted scattering will give some hints at new physical opportunities that can be expected in 
collisions of high-energy particles carrying OAM, while the calculation of a twisted particle 
decay clarifies various technical details involved in passage from plane waves to twisted states. 

The OAM and spin are two forms of angular momentum, and the problem of gauge- 
invariant separation of these two objects has a long history. It is still being debated both 
in the optics community, see for example [10], and in the HEP community, especially in 
the context of the notorious proton spin puzzle, [TTj . For the problems we consider here it is 
sufficient to note that all experimental situations which seem to be realizable are well described 
by the paraxial approximation. It is known that in the paraxial approximation spin can be 
well separated from the ^-component of OAM of light, [12]. The same applies to fermions as 
well. Therefore, incorporation of both spin and OAM degrees of freedom, leading to non-trivial 
polarization fields, does not seem to pose any problem in the paraxial approximation. 

However in the present paper we would like to focus specifically on the OAM component 
and to understand what new physical opportunities are offered by the non-trivial spatial 
dependence rather than unusual polarization states. Therefore we limit ourselves to scattering 
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of scalar particles only. In this simple case all the non-zero angular momentum is definitely 
due to the orbital part. Additional features arising from the polarization parameter fields will 
be considered separately. 

The paper is organized as follows. In Section |2] we introduce scalar twisted states and 
describe some of their properties. In Sections [3] and H] we derive expressions for the cross 
section in the single-twisted and double-twisted cases, respectively. Section |5] gives a thorough 
discussion of kinematical features arising in the case when the final state contains twisted 
particles. In Section [6] we discuss the results obtained and draw our conclusions. In two 
Appendices we derive some technical results used in the paper. 



2 Describing twisted states 
2.1 Spatial distribution 

As mentioned in the introduction, we focus in this paper on twisted scalar particles with mass 
M. In their description we follow essentially [HIE]. 

We represent a state with a non-zero OAM with a Bessel beam-type twisted state. This 
is a solution of the wave equation in the cylindric coordinates with a definite energy u and a 
longitudinal momentum k z along a fixed axis z, a definite modulus of the transverse momentum 
|k| (all transverse momenta will be written in bold) and a definite ^-projection of OAM. If 
the plane wave state |PVK(k)) is 

\PW{k)) = e - iLUt+ik * z ■ e ikr , (1) 
then a twisted scalar state \k, m) is defined as the following superposition of plane waves: 

-lUlt+lk z Z f d k flrW kr r, flr\ — ( n\ m a im 4>k ^ /oZ^(W ^) 



\ K) m) = e -"*+»*.* / QKm ( k ) e ^ , Gftm ( k ) = (-t)*V"0»V27r Vl 'j= ' . (2) 



In the coordinate space, 



e im<f> r 



\K,m) = e- Mk * z -iP Km (r), Vw(r) = -=V^J m («r) . (3) 

V Z7T 



Here, following [5] we call k the conical momentum spread, m is the ^-projection of OAM, 
and the dispersion relation is k^k^ = u> 2 — k\ — k 2 = M 2 . We note in passing that the average 
values of the four-momentum carried by a twisted state is 

(P) = (w, 0, k z ) , (4) 

so that (k' 1 )(k l j) = M 2 + k 2 , which is larger than the true mass of the particle squared. 
The transverse spatial distribution is normalized according to 

- t,^ Jr^U^r) = *^(. - V) . (5) 
The plane wave can be recovered from the twisted states as follows: 

\PW(k=0)) = limW— k,0), (6) 
«->o V k 

|Piy(k)) = W— i m e~ im(pk \K,m) , «=|k ± |. (7) 

m=—oo 
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If needed, these two cases can be written as a single expression: 



\PW(k)) = 




From these expressions one sees that the twisted states with different m and k represent 
nothing but another basis for the transverse wave functions. 



dn PW = tiR 2 L z /o 2 . (9) 



2.2 Density of states 

When calculating cross sections and decay rates, we need to integrate the transition proba- 
bility over the phase space of the final particles. When calculating the density of states, we 
consider a large but finite volume and count how many mutually orthogonal states with pre- 
scribed boundary conditions can be squeezed inside. In the present case due to the cylindrical 
symmetry of the problem, we choose a cylinder of a large radius R and a length L z . In the 
case of plane waves we have 

dk z d 2 \t 

w? 

The full number of states with transverse momenta up to kq and longitudinal momenta \k z \ < 
k z0 is k z0 L z ■ R 2 nil 'Air. 

To count the number of twisted states \k, m) in the same volume, we specify the boundary 
condition, e.g. ip Km (r = R) = 0, which makes k discrete such that KiR is the i-th root of the 
Bessel function J m . Then we note that the position of the first root of the Bessel function 
J m {x) is always at x > m, and as m grows x — > m. For a given k, the maximal m for which 
the wave can still be contained inside the cylindrical volume is m max = kR, which has a very 
natural quasiclassical interpretation. 

If m is small and not growing with R, then one can use the well-known asymptotic form 
of the Bessel functions to count the number of states: 

RdnL z dk z Am 
dn tw = — . (10) 

Here, Am is written instead of just 1 to signal the presence of a discrete running parameter 
m. 

If m is not restricted to small values, this asymptotic form of J m (x) cannot be used since 
it requires m? <x. Instead, the so-called approximation by tangents can be used, which gives 
the following density of states: 

-dn Am L 7 dk ? 



dn tw = v 7 m max ~ ™ 2 TT^ ■ (H) 

K 71 Z7T 

In the limit m m max = k,R this expression reproduced fflUj) . Alternatively, one can calculate 
the radial part of the density of states via the adiabatic invariant as suggested in |6|. The 
number of radial excitations n r for a fixed m is 



f R k r (r)dr , . . / „ w? , , 

n r = ^L^ ; kr (r) = JK 2 - — . (12) 
Jm/K n Mr 2 
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The density of states is then given by 



dn r = = \/ m 2 max - m 2< ^- . (13) 

One important remark is in order. Effectively, switching from the plane wave to twisted 
state basis for the final particles implies replacement 



, 9 , / m 2 , Am . . 

d 2 k^4Wl —ndn . (14) 

V m max m max 

Note that the contribution of each "partial wave" with a fixed m vanishes in the infinite 
volume limit as 1/R. However, the number of partial waves grows oc R, and in order to get 
a non-vanishing result for a physical observable, one must integrate over the full available m 
interval up to m max . This remains true even if the transverse momenta stay small, and it is 
related to the fact that the plane wave contains contributions from all impact parameters with 
respect to any axis non-collinear to its propagation direction. 

Another expression one needs for the probability calculations is the normalization constants 
for the one-particle states. A usual plane wave one-particle state is normalized to 2E ■ V; to 
renormalize it to one particle per the entire volume, the plane wave should be multiplied by 
N PW) with 

N 2 PW = ^ V , V = nR 2 L z . (15) 
For a twisted state the corresponding normalization factor N tw is 

9 1 7TK, , . 

i\& = —-=_=-, (16) 

which in the small-m case simplifies to 

j\r 2 «- — (17) 

also derived in [B] . Note however that even in the general case the product of the normalization 
constant squared and the density of states for each final twisted particle is simplified as 

NlMt.. 



"tw u '"'tw 



2E ■ 2vr 



2.3 Flux and the cross section 

The definitions of the flux factor and the cross section have to be reevaluated when a collision of 
non-plane wave states is considered, which involves subtle issues described in [6] . By definition, 
the cross section is the transition probability per unit time divided by flux. In the plane wave 
case, when the four-momenta of colliding particles are fixed, both the flux and probability 
are constant across any chosen plane, so the proportionality between them holds locally. For 
a head-on collision one gets jpw = (|fi| + 1^2 D/V', which together with the energies of the 
incoming particles and one volume factor combine to the familiar Lorentz invariant expression 

I PW = V(pk) 2 ~P 2 k 2 . (19) 
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This formula can of course be used in any frame, including cases when the collision is not 
head-on. 

In the single-twisted case, when a twisted state collides with a plane wave, both the flux 
and the transition probability are not constant but change across the transverse plane. As 
noted in [6], one therefore needs to redefine the notion of the cross section to adapt to this 
situation. One introduces the averaged cross section defined as the transition probability 
integrated over all r divided by the flux again integrated over all r. While the definition of the 
former quantity is clear (this is what we calculate in the following two Sections), the proper 
definition of the integrated flux is more intricate and apparently not unique. A definition of 
the flux should however be correlated with a definition of the (generalized) luminosity, so that 
the observable event rate remains uniquely defined. 

In [6] the forward electron-photon collision was considered (the 3-momentum of the plane 
wave electron was directed exactly along the axis z), and the following procedure was sug- 
gested: the total flux is just sum of the z-components of the two fluxes (note that we actually 
assume the absolute values of the fluxes): 

«>=« + <C> = ^^, (20) 

where tano^ = n/k z . We think that a more appropriate definition of flux should take into 
account not only the ^-components of the individual fluxes but also the relative lateral motion 
of the two waves. Indeed, the sole purpose of calculating the flux factor is, classically speaking, 
to derive the volume swept by one particle in the "gas" of opposing particles and find how 
many "attempts" at collision are made per unit time. Therefore, we propose the following 
general definition of the integrated flux factor I tw for the single-twisted case: 



' tw 



2tt 



fc W(k,p). (21) 



Note that it is well defined for any transverse momentum p of the opposing plane wave. In 
the specific case considered above this definition gives (j) = (1 + v cosak)/V, which differs 
from (|20|) . This definition can be also generalized to the double- twisted case: 



L2tw 



2tt 2tt 



W(k,p). (22) 



3 Single-twisted cross section 

We start with a usual 2 — > n collision in which both incoming particles are described by plane 
waves with definite four-momenta k and p. The final system X is also treated as a collection of 
plane waves with the total momentum px- The invariant amplitude of this process is denoted 
by .M(k, p), where the transverse momenta of the initial particles are indicated explicitly. 

The cross section of this process is calculated according to the standard rules. When 
squaring the scattering matrix element 

S PW = i(2n)W A \k + p-p x )- M{k, p) , (23) 

we re-interpret the square of the delta-function as 

[5 (4) (k + p - p x )] 2 = (k + p-p x ). . (24) 
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Using the plane-wave normalization factors ( 11 5 p for all the initial and final particles, we get 
the cross section 

M , P) = (2 ' rm - EMP ;- P " mk + P - PX) \M(k, p)|- . dT, . (25) 



PW 



For future convenience the delta-function is explicitly broken into the longitudinal and trans- 
verse parts. As usual, we can extract from the final phase space integration measure dY x the 
integral over the total transverse momentum px, dT = d 2 pxdT' x , and write the cross section 
as 

a fi \ {2n) A 5{Ei- E f )5{p zi -p zf ) 2 , 

da PW (k,p) = — 1 ^|.M(k,p)| -dT x . (26) 

4i pw 

Now we recalculate the cross section for the case when the first particle is in the twisted 
state | «, m). We apply prescription (j2J) to the ^-matrix element: 

/d 2 k 
— —a Km (k)S PW (k, p) , (27) 

as originally suggested in [5]. The plane- wave S- matrix ( 123|) contains a delta- function of 
transverse momenta, which makes it possible to simplify the square of the twisted S'-matrix 
element as 



f crk crk 

Stw? = j j^j^a Km (k)a* Km (k')S PW (k,p)S PW (k',p) 

oc J ( 2^ 2 0j2 a U*-K m (k')8W(k + p - p x )5^(k' + p - p x )-M(k, p).M*(k', p) 
l (k)< m (k)^(k + p - px)|M(k, p)| 2 . (28) 



(2vr) 4 

The square of a Km (k) contains a radial delta-function squared, which is reinterpreted as 

[SU - |k|)] 2 = 5(k - |k|) ■ 5(0) — >• 5(ft - |k|)- . (29) 

7T 

This prescription comes from the observation that at large but finite R and at k = \k\ the 
radial delta-function is regularized as 

5(0)= rdr[J m (Kr)} 2 -)■/ rdr[J m (Kr)f w - . (30) 



7T 



see [3 E]. Therefore, with all the normalization factors ( |T5l) and ( 1T7|) the single- twisted cross 
section takes form 

f d 2 k I PW (k,p) 5(k - \k\) f d(j) k I PW (k,p) 
d °tw= / ~z ; da PW (k,p) = / da PW (k,p) , (31) 

J 27T l tw K J Z7T J tw 

see Section 12.31 for the definition of the fluxes. Note that in the paraxial approximation one 
can replace the ratio of the fluxes by the unity. 

A couple of remarks concerning this result are in order. In the usual case of a plane wave 
collision, the total initial and, therefore, final momenta are fixed. This is highlighted by the 
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presence of the transverse delta-function in (j25p and by implicit correlations among transverse 
momenta of the final particles inside dT' x in fl26|) . In the process involving a twisted particle, 
the total final momentum is not fixed, making the final particles less correlated. Although 
all the plane waves |PW(k)) which constitute the initial twisted state \k, m) are summed up 
coherently at the amplitude level, the final states with different momenta do not interfere, 
and this breaks the coherence. As a result, the twisted particle cross section is expressed via 
an incoherent superposition of the cross sections induced by each initial plane wave. This is 
precisely what (|3T|) displays. 

The previous paragraph contains one additional subtlety. When saying that the final states 
with different momenta do not interfere, we implicitly assume that the final state is detected 
by a usual detector, which measures the linear momentum but not the OAM. On the contrary, 
if the final particles were detected by a hypothetical "coherent detector" sensitive to a coherent 
superposition of final states with distinct momenta, or if the production of particles in the 
reaction we consider is followed by another process which is OAM-selective, then the coherence 
would be restored. Thus, the coherence is not actually destroyed in the scattering process, 
but remains hidden. 

The second remark concerns the angular region contributing to the integral (l3Tj) . In the 
fully differential case, that is when we fix the momenta of all the final particles, the transverse 
delta-function inside dapw assures that there is only one value of (f) that contributes to the 
integral. At this level of consideration, representing the cross section as an angular integral 
might look somewhat misleading. However if this delta-function is killed by the integration over 
all allowed px, or equivalently by integrating over one of the final particles, then the integral 
receives contributions from all angles Note that the angular integral representation is 
also justified even in the fully differential case if the detector resolution of the final particles' 
momenta is worse than k. 

Let us also show a slightly different derivation of ( 1311) . We first explicitly perform the 
integration in (|2"T1) . which is effectively killed by the transverse part of the delta-function, 
keeping the scattering amplitude essentially unchanged: 

/ ^ a ^(k)S^(k + P-Px)-M(k, P ) = y^e im *« 6 -^fi ■ P ) , (32) 

where q = px — P- Note that after the integration the transverse momentum k cannot appear 
in the amplitude any more. However, since the other particles are plane waves with well 
defined momenta, the vector q with modulus q and azimuthal angle <p q is also well defined 
and can be used instead of k everywhere. Squaring (1321) . we again encounter the square of the 
radial delta-function which is reinterpreted as in (1291) . The cross section then becomes 



. (2it)%E t -E f )5( Pzt - Vzf )\M^p)\ 2 6(k - q) 
do tw = — al x ■ — ^ • (33) 

The integral over the overall transverse momentum present in dTx kills the radial delta- 
function: 

2 8{K-q) f 2 S(K-q) f d<j) q 



2ttk / ZHK / 27T 



and one recovers the result fl3Tl) . 

One can make several observations concerning fl3"T|) : 
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• The cross section is m-independent. 

• The initial twisted particles effectively perform the angular averaging of the plane wave 
cross-section. 

• There is no smallness associated with non-zero m. 

• There is no small factor associated with small k. 

Let us now consider the case when the twisted particle is not in an eigenstate of the OAM 
operator, but in a superposition of states \k, m) with equal k but different m. For example, 
consider the twisted state of the form a\n, m) + a'\n, m'), with \a\ 2 + \a'\ 2 = 1. Repeating the 
same calculation we encounter an interference term in the cross section: 

a Am f d-4>k Ipw(k, p) , . , . , . 

da tw = / n r cos(Am0 fc + a)dap W (k,p) , (35) 

J 2n Uw 

where Am = m — m! and a is the relative phase between the two complex coefficients a and 
a'. The cross section for such an initial state takes the following form: 

da = da tw + 2\aa'\ da^ 1 . (36) 

Results (EH) and ( 1351) mean that if the initial particle can be prepared in a twisted state 
with an adjustable superposition of different m, a Fourier analysis of the cross section with 
respect to the initial azimuthal angle can be performed. In principle, the same analysis can 
be done with plane waves, but it would require making several experiments with different 
angles of the initial particle k and then extracting the Fourier components via the partial 
wave analysis. From the experimental view, it is likely that systematics of the two schemes 
can be different, which makes them complementary to each other. 



4 Double-twisted cross section 

Let us now consider collision of two initial twisted particles. In this work we do not aim at a 
systematic study of this case, but rather outline some new features can be expected in such 
circumstances. Therefore we consider the simplest set-up, in which the two colliding particles 
are described by twisted states \k, m) and \r), n) defined with respect to the same quantization 
axis z: 

km) = f ^^(kJlPWiCk)) , \ V ,n) = J ^- 2 a m (p)\PW 2 (p)) , (37) 

with the same functional form of the projectors a as before. Here, the subscripts 1 and 2 refer 
to the first and second colliding particles. The "double-twisted" version of (!2Tj) is 

72 \r Al 



/d k dp 
7^—^ 77TT5 a Krn(k)a vn (p)Sp W (k, p) , (3£ 
(27r) 2 (2tc) 2 



and its square is proportional to 

d 2 kd 2 pd 2 k'd 2 p' 



( , )t) s a™(k)a,„(p)o* ra (k>yp') 
x <^(k + p - p x )S {2 \k' + p' - p x )M(k, p)Af(k', p') . (39) 
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Let us consider the kinematical restrictions imposed by the delta-functions entering this ex- 
pression and compare them with the result k = k' found in the previous Section. Here, too, 
the moduli of the momenta are fixed and pairwise equal: |k| = k — |k'| and |p| = 77 = |p'|. 
Besides, the two pairs of momenta sum up to a well-defined px- For each pair there are two 
possibilities satisfying these conditions, shown in Fig. HJ which are mirror reflections of each 
other with respect to the direction of px- Therefore, the integral ( |39|) receives contributions 
from two kinematical configurations: 

direct: k' = k , p' = p , 
reflected: k' = k* = -k + 2(kn x )n x , p' = p* = -p + 2(pn x )n x , (40) 

with nx = px/|px|- As we will see below, in contrast to the "single-twisted" case, here the 
existence of two configuration for any given px 7^ leads to a residual coherence between 
different plane wave components in the twisted state. 




Figure 1: Two kinematical configurations of the transverse momenta k and p of fixed absolute 
values that sum up to the vector p X - 

Our goal now is to see whether the double-twisted cross section can be expressed as an 
angle-averaged plane wave cross section similarly to ( I3~TT) . To this end let us first study a 
generic expression 

J = Jd<p k d<p p 5^(k + p-px)-f(k,p), (41) 

where |k| = k and |p| = rj are fixed and /(k, p) is a continuous function of the incoming 
particles' momenta. Due to the delta-function, the integral J receives contributions only from 
two points, shown in Fig. [IJ At these points the azimuthal angles 0^ an d 4> p take specific 
values: 

4>k = 4>x ± S k , <f)p = 4> x =F S p , (42) 
where (px is the azimuthal angle of px and 

f p 2 x + K 2 -r] 2 \ f Px~ K 2 + r] 2 \ 

dk = arccos ; : , dp = arccos ; ; (43) 

V 2|px|/« J V 2|p x |r? J 

are functions of the absolute values of the momenta. Let us denote the values of /(k, p) at 
these two points as /+ and /_, respectively. In Appendix [A] we derive the following result 

J = f - ± ^~ = J dfa d^ 2 \k + p - px) ■ f - ± ^ , (44) 
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where A is the area of the triangle of sides |px|> K i V- ^ allows us to rewrite the square of the 
integral J as 



\j\ 2 = y# fc # P ^(k+p- PX )/(k,p)- 



ft+f- 



2A 

= ±J ^ fc ^^)(k + p-p x )[|/ + | 2 + |/_| 2 + 2Re(/ + /:; 

= d<M<M (2) (k + p- P x) 0/(k,p)| 2 + Re/(k,p)r(k*,p*)] , (45) 

where k* and p* are given by (j40l) . 

These results can be directly applied to the integral (EHD if we note that it has the form of 
| J | 2 with the function 

/(k - p) =(2^ e ™'- M(k - p )' 
The integral (159"j) then takes form 

1 



J d0 fe #p5 (2) (k + p-p x ; 



(27r) 6 sin(4 + 5 p ) 

x {|7W(k,p)| 2 + Re [e 2im( ^-^ )+2in( ^- fe) >l(k,p)7W*(k*,p*)]} . (46) 

Bringing together all the normalization coefficients, we finally arrive at the following repre- 
sentation for the double-twisted cross section: 

da 2tw = q — r-h —ft / rf0 fc t% Ipw ^ p ) [da P vy(k,p) +dtr'(k,p)] , (47) 

57T SUl^Ojfe + OpJ J l2tw 

where the flux I 2tw is given by (1221) and 

da'(k,p) = ( 27r )^(^- ^/)^-^/) Re [ e 2m (^-fe)+2,n(0 p -^) A/1(k)p) _ M * (k * )P * ) ] . ^ 

4:IpW 

(48) 

We see that the "direct" contribution yields the standard cross section, while the "reflected" 
contribution gives rise to the novel quantity <x'(k, p). This quantity describes the auto- 
correlation of the amplitude and is absent in the plane wave case. 

Let us also see how ( 14?]) simplifies in the case when .M(k, p) = .M(k*,p*). The novel 
quantity reduces to the usual cross section, da' = cos[2(m<5fc — n8 p )]dapw, and we obtain 

cos 2 {m5 k -n5 p ) f I PW (k,p) 

da 2tw = — — — / d(p k d<p p da PW [k, p) . (49) 

47rsm(() fe + 5 P ) J htw 

It is only in this case that the double-twisted cross section is expressed via the angular integral 
of the plane wave cross section. 

As in the single-twisted case, we note that if all the final momenta were fixed and if 
the detector had an infinitely good momentum resolution, there would be just two points 
contributing to the angular integral (147p . If at least one of these conditions is broken, the 
integrand extends over the full integration range. 
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Similarly to the single-twisted case, the double-twisted cross section stays finite even when 
the conical momentum spreads k and rj are very small. However in contrast to the single- 
twisted case, the cross section now depends also on the orbital angular momenta m and n. 
This dependence comes solely from the correlation term in ( 1471) . If the cross section is averaged 
with a sufficiently smooth function of the initial or final momenta, then in the limit of large m 
and n the correlation term gets suppressed by rapid oscillations, and the double-twisted cross 
section stays approximately m and n independent. 

Certainly, one can also consider double-twisted cross sections with twisted states in su- 
perposition of different m and n. We leave a detailed study of this situation for a future 
work. 

Finally we reiterate the point that for this particular version of double-twisted scattering 
we chose the simplest possible set-up, in which both twisted states are defined with respect to 
the same z axis. One can also study what would happen if two different axes, either parallel 
or not, were used. Answering this question requires a more elaborate formalism. 

5 Two-particle decay of a twisted scalar 
5.1 Preliminary remarks 

Previous two Sections were devoted to the cases when the twisted particles appeared only 
in the initial state, while the final state X was assumed to be describable with the plane 
waves. Let us now consider the case when at least one of the final particles is also twisted. 
This is exactly the situation that arises in the original suggestion of E] to use Compton 
backscattering of the twisted optical photons to produce final high-energy OAM photons. The 
key question here is to what extent the twisted parameters of the initial state \k, m) determine 
the parameters of the final twisted state \K',m'). Here again both twisted states are defined 
with respect to the same z-axis. In the case of strictly backward scattering the conservation 
of the orbital momentum ensures that k' — k, m' — m, [5]. In this Section we study how this 
result changes for the non-forward scattering. Note that we use the term "forward" for any 
process with zero transverse momentum transfer, i.e. both for strictly forward and strictly 
backward kinematics. 

In order to focus on the final state kinematics and to avoid possible complications coming 
from a non-trivial matrix element, we would like to answer the above question in the context 
of the simplest possible problem: the decay of a twisted scalar particle with mass M into a 
pair of massless distinguishable particles due to the cubic interaction g ■ $0!0 2 . The momenta 
of the initial and final particles are p, k\ and respectively. We will calculate the decay 
width in the center of mass frame defined by p z = 0. This is not the true rest frame because 
due to the transverse motion a twisted particle is never at rest. To make the presentation 
more pedagogical, we will first calculate the decay rate when both particles in the final state 
are plane waves, then for the plane wave plus twisted final state, and finally for the case when 
both final particles are twisted. Although the total decay width must be the same in all these 
cases, the differential decay rates will be rather different. 
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5.2 Two plane waves 

Again, let us first recall the standard calculation for the case when all the particles including 
the initial one are plane waves. The S'-matrix is given as usual by S = i{2t^) a 8^ (p — ki — k 2 )-g. 
With the plane wave normalization coefficients for all the particles, the differential decay rate 
for a particle at rest is 

dT = MV* W (P - h ~ k ^ VT ■ {N l w f ■ dnpw (h)dn PW (h) 

- 92 6{M -^-^k u (50) 



(2tt) 2 

so that the total width is 

Now we repeat this calculation for the initial twisted state \k, m), while keeping the plane 
wave basis the final particles. The S'-matrix is 

/72l 
^a Km (k)5( 2 )(k-k!-k 2 ) (52) 

= i{2*) A g6{E -ujx- u 2 )5(k lz + k 2z )^^e im ^ K ^ 12) , (53) 



where k\ 2 = \/k 2 + k| + 2|ki||k 2 | cos(0i — 2 ) an d <pi 2 is the angle of the 2D vector k 12 w.r.t. 
some axis x. Similaly to the single scattering cross section, this phase factor is inessential and 
m disappears in the decay rate. 

The square of 5(n — k\ 2 ) is treated as in fl29l) . and with the appropriate normalization 
factors taken into account, the decay rate has the form 

(2tt)V 5(K-k 12 )R 1 vr Vd 3 hVd 3 k 2 

dT = — —5(E -u)i-uj 2 )6(k lz + k 2z )TL 



ZEuj^-T v ' v z z> z k vr V 2 RL Z (2tt) 3 (2tt) 3 

g 2 8(k - k 12 ) 5(E - u x - u 2 ) „ , 2l _ j2 . 



(2tt) 3 k SEu^ 
For the transverse integral we write 



dk z d 2 k x d 2 k 2 . (54) 



J d<j> 2 ^ ^ 12) = 2 J d<f> 2 6 [k 2 - k\ - kj - 2|ki||k 2 | cos(0! - <f> 2 )] = i , (55) 

where A is the area of the triangle with sides k, |kx| and |k 2 |, see Appendix [A1 As usual, the 
energy delta function can be killed by the k z integration 

uj\uj 2 Ek* z 
where 

1 
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The decay rate becomes 



A • (58) 



The integration region over |ki| and |k 2 | is defined by the requirement that, in addition to 
(j8ip . the longitudinal momentum k* z is well-defined, which cuts the rectangular shape shown 
in Fig. [2] in Appendix |A] It is conveniently described with variables 

|ki| — | k2 1 |kx| + |k 2 | __ E 

x , z , x G [ 1, 1J , z G [1, z max ^ , z max = 1 • (59) 



hi 



In these variables 



a = ^v / ^-i)(i-- 2 ), k : = jJ(i-^-)(i-^), (60) 



max / \ max 



and the decay rate takes form 



dT = -4 (^-x^dx 

[(4 a ,-^)(^-^)(^-l)(l-^ 2 )] 1/2 
This integral can be taken exactly, and it gives the decay width of the twisted scalar particle 

(62) 



16ttE 



which is a very natural result. In the limit k — y 0, we recover the plane wave decay width 

(ED). 

Of course, one could also arrive at (1621) just by using our previous results concerning the 
single-twisted cross section modified to the case of one initial particle. However the detailed 
derivation given here is needed to understand what changes when the final state includes 
twisted particles. 

5.3 Twisted state plus plane wave 

Let us now describe the final state as twisted state plus a plane wave: 

|«, m) -> |«i, m x ) + |fW(k 2 )> . (63) 

Since the full decay width cannot depend on the basis we choose for the final particles, we 
must recover the same result fl62|) in this basis. In addition to that, we also want to know how 
the final twisted state parameters k% and mi are related to the initial state parameters k and 
m. 

The S'-matrix is now 

S = J (2^)2 (2n) 2 a * Kimi ( kl ) flKm ( k ) Spw ■ ( 64 ) 

It can be written as 

S = i(2Txfg5(E -u x - u 2 )5{k zX + k z2 ) ■ X m>mi («;, «i,k 2 ) . (65) 
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where the master integral I m mi (/t, k 1; k 2 ) is defined as 

72 U j2 



/d k cfk 
7^72^)2<im 1 ( k l) a «m(k)5 (2) (k - ki - k 2 ) . (66) 

Let us first calculate this integral in the strictly forward case k 2 = 0: 

x„(^ 1; o) = / d 2 kd 2 k ie ^-^^ ^^^4^ ^)(k-k0 

= y # #1 e iTn *- imi * ■ 25(k 2 - k 2 )5(</> - 00 

- Ki)S m>mi , (67) 



(2vr 

which was first obtained in [5] . This result implies that the twisted state quantum number are 
transferred from the initial to the final twisted particle without any change. The differential 
decay rate is 

" r = 37^ o m>rni d{K - Ki) ■ dKidk lz d k 2 

(2ti) 2 8EU1OJ2 



g 2 5{E-u l -u 2 ) ,3, 9 d k - 



d 6 ki= » A ' , (68) 



(2tt) 2 8E0J1OJ2 (2ir) 2 4E 2 k 

and we stress that this result is applicable only at k 2 = 0. 

In the general non-forward case the master integral ( 1661) can be rewritten as 



•mi-m r 

l m , mi (K } K U k 2 ) = y ##ie imHmi ^ -^(k-ki-ka). (69) 

There are two ways to look at this integral. First, using the results of Appendix [A] we obtain 
W,*,*) = ^VSfr«— » "' ImA ~ f " "" W , (70) 
where A is the same as in flHBI) with k 2 replaced by Additionally, we can also rewrite 

<5 (2) (k - ki - k 2 ) = — / d 2 r e irk - irk i-' rk 2 (71) 

and represent the master integral as an integral over a triple product of Bessel functions which 
is useful in certain circumstances: 

jmi—m e i(m—mi)$2 roo 

2m,mi(K,Ki,k 2 ) = — — ^[kk[ ■ I rdr J m (/cr) J mi («ir) J m - mi (|k 2 |r) , (72) 

where 2 is the azimuthal angle of k 2 . Note that although we are considering the case with 
two twisted particles (one in the initial and one in the final state), an integral over three Bessel 
function arises automatically. Comparison of (!72|) with (1701) gives the result for the integral of 
the triple Bessel function product. See also [13] for some mathematics involved in evaluation 
of this and similar integrals. 
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Let us also check the k 2 — > limit of the result ( 17D|) . When |k 2 | <k, the distribution over 
Ki spans from k — |k 2 | to k + |k 2 |. The angle 5% — >■ 0, while 5 2 can still be arbitrary. However, 
in the k 2 — > limit only m — mi = term survives due to J m _ mi (|k 2 |r), so that the cosine in 
( !?0l) approaches unity. The analysis of A shows that 



lim = 2tt5(k - «i) , (73) 

|k 2 |->o A v 1 



so that one indeed recovers the strictly forward result for the master integral given in (J6j 
Note that although we wrote simply |k 2 | — > 0, we imply a very specific combination of this 
limit with the limit R — > 0, an issue to be discussed in Section 16.21 

We continue the calculation of the decay rate in the non-forward case. After integration 
over k z , the decay rate can be written as 

dT = 47TV ^ " %n, mi (K, «*, k 2 )| 2 • ^ d 2 k 2 . (74) 

Note that this decay rate is differential not only in K\ and k 2 but also in the discrete variable 
mi, the full decay width includes integrals over momenta and a summation over all possible 
mi's. 

A close inspection shows that the immediate integration over K\ or k 2 cannot be done due 
to singularities of |Z| 2 along the boundaries of the kinematically allowed region. In contrast 
to the plane waves in the final state, the denominator now contains A 2 instead of just A. 
Therefore, in terms of variables x and z one encounters singularities of the form 

dx . [ Zmax dz 

and 



x 2 /i z 2 



Clearly, this is an artefact of the infinite radial integration range in (J72J). If instead we take R 
to be large but finite, we expect that a trick similar to (1291) should be at work, namely that 
after regularization |X| 2 would yield R times a less singular function. 

This trick does not seem to work for each mi separately. However, as we prove in Ap- 
pendix (Bj it works for |I| 2 summed over all possible mi. In the limit R — > oo we obtain 

^ |X mimi ( K)fi;i ,k 2 )| 2 = ^^^-, (75) 



m\=— oo 



with the same A as before. The regularization parameter R then disappears from the result, 
and the decay rate reads 

,-p 9 2 1 KidKi |k 2 |<2|k 2 | 

dT = ^T 2 4E 2 X A • (76) 

Comparing (JTB"]) with the previous results (|68|) and fl58l) leads us to two conclusions. 

• The transition from the strictly forward to the non- forward cross section/decay rate 
consists in replacement 

_1 (77) 
m 2vrA v ; 
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• The result (I75|) for a twisted particle in the final state coincides with the result (I58p for 
the case when both final particles are plane waves. 

Although these conclusions were drawn for the specific process we consider, it represents a 
universal kinematical feature characteristic to all processes involving non-forward scattering 
of twisted states. 



5.4 Two twisted states 

For completeness, let us also recalculate the decay rate in the basis when both final particles 
are described by twisted states: and \k 2 , m 2 ). We remind that all twisted states are 

defined with respect to the same common z axis. It turns out that this calculation closely 
follows the case of twisted state plus plane wave just considered. This is not surprising because 
the appearance of the triple Bessel integral highlights the fact that when two particles (one 
in the initial and one in the final state) are twisted, the third one is automatically projected 
from the plane wave onto an appropriately defined twisted state as well. 
The S-matrix takes the form 

— ui — L0 2 )S(k z i + k z2 ) ■ <5 m , mi+m2 v /KK i K 2 J rdrJ m (Kr)J mi (Kir)J m2 (K 2 r) , 

(78) 

and we again encounter the triple Bessel-function integral. The decay rate is written as 

47r 2 4£ 2 A;*' A 1 } 

This expression is identical to (ITS]) up to the obvious replacement | k2 1 — > k 2 ; its integration 
over all conical momenta spreads K\, k 2 gives again ( 1621) . 



6 Discussion 



6.1 Potential uses of the twisted cross sections 

In this paper we represented the cross sections for the single-twisted and double-twisted pro- 
cesses via the angular integrals f )3T|) and ( 147]) . These are rather unconventional quantities, as 
they involve averaging of the plane wave cross sections over the azimuthal angles of the initial 
particles at fixed final momenta. A further study is needed to see if they can be related to 
the more conventional cross sections, in which the initial momenta are fixed but final state 
azimuthal angles are integrated out. 

Using twisted states in superpositions of different orbital angular momenta, one can per- 
form, in principle, a Fourier-analysis of the differential cross section as the function of the initial 
azimuthal angles at fixed final momenta, see (1351) . This might be a useful complementary tool 
for the study of various azimuthal asymmetries. 

Perhaps, the most intriguing feature we have found arises in the double-twisted cross 
section. We showed that this cross section is sensitive not only to the plane wave cross 
section, but also to a quantity describing the autocorrelation of the plane wave amplitude. 
Such quantity is absent in the conventional plane wave scattering. 
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Note that it is not the first time the amplitude autocorrelation appears in particle physics. 
In [H], where finite beam size effects were analysed, the product of amplitudes at different 
momentum space points was also involved. The essential novelty of the present case is the 
extra angular factors in ( )48|) . which offer a more detailed probe of the amplitude. 

Finally, let us stress once again that experimental study of these features does not require 
large conical momentum spread k and is feasible even for k in the eV range, which was already 
realized for energetic electrons in [Hi |9] and which was proposed in [5] for high-energy twisted 
photons. It is the non-trivial angular dependence that is at work. However we stress that in 
order to probe the autocorrelation function of the amplitude with the momentum offset q, one 
would need conical momentum spread k ~ |q|. 

It is natural to ask whether the new degree of freedom offered by twisted particles, and in 
particular photons, can be used to gain more insight into the structure of hadrons. Indeed, 
the structure of proton remains one of the hottest topics in QCD. One of the questions that 
might be particularly relevant for the present study is the contribution of the quark/gluon 
orbital momentum to the proton spin, see e.g. reviews [UJ. It is, therefore, very interesting 
to know whether the OAM of the initial photon can be transferred to partons, allowing us to 
probe the partonic composition of the proton in a novel way. 

Without going into detail, let us make some preliminary observations that can be inferred 
from the results of the present work. 

Let us first note that the definition of a twisted state ([2]) involves an angular functional: 
it picks up a plane wave quantity as a function of the initial azimuthal angle 0, weights it 
with e tm ^ and integrates over all angles. This might give an impression that the twisted state 
is a perfect Fourier analyser at level of amplitude. This impression is wrong because this 
functional acts not on the amplitudes itself, but on the S- matrix element, which includes a 
transverse delta-function. Eq. fl32|) shows that this functional partially kills the delta-function, 
while the amplitude remains unchanged. The orbital angular momentum is just transferred 
to the overall motion of the final system instead of exciting internal degrees of freedom. This 
somewhat pessimistic conclusion is supported by the result of this paper that the single-twisted 
cross section is just an azimuthal integral of the plane wave cross section. 

The situation might change if both the initial photon and proton were twisted. As we 
showed above, in this case a residual coherence between two kinematical configuration survives, 
which leads to an additional term in the cross section proportional to the autocorrelation 
function of the amplitude. This term has a non-trivial azimuthal structure and might represent 
a new probe of the proton. This issue certainly deserves a dedicated study. 

6.2 Properties of the final twisted particle in a non-forward scat- 
tering 

Let us also discuss in detail what the results of Section 15.31 tell us about the values of Hi and 
rrii. The differential decay rate f !76p shows that at large |k 2 | ^> n the conical momentum 
spread K\ is limited to the interval from |k 2 | — k to |k 2 | + n with an inverse square root 
singularity at the endpoints. This singularity is integrable, so that the entire interval more or 
less homogeneously contributes to the integral. Since |k 2 | can be as high as E, the total decay 
width is therefore dominated by large K\ ^> k. In a more complicated process, the decay rate 
or the cross section will include the amplitude squared which can serve as a cut-off function. 
For example, in the Compton scattering one expects that K\ up to ~ m e will contribute to 
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the cross section. 

Although we found a result for the decay rate summed over all m 1; we can trace the main 
mi-region from the intermediate formulas. The result is that essentially all mi from minus 
to plus infinity are important for the decay rate, which is in a strong contrast to the strictly 
forward result m — mi. 

Indeed, the distribution over final m\ can be seen in (I74p . where the master integral I m ,mi 
is given by (ITU]) . For generic transverse momenta, growth of m\ leads to oscillations of the 
cosine function with constant amplitude. Thus, when averaging over the entire Ki interval, 
one can approximate cosine squared by 1/2, and the dependence on mi drops off. This result 
holds for any non-zero transverse momentum transfer |k 2 | —>•(). 

Since the forward and non- forward mi-distributions are so dramatically different, a natural 
question arises whether there is a continuous transition from the non-forward to the forward 
scattering. The answer to this question involves an accurate treatment of two limits: | lc2 1 — > 
and R — > oo. Let us keep R large but finite, and set |k 2 | — > 0; then the transition is 
smooth. Looking at the triple-Bessel representation of the master integral f lT2"j) with the upper 
limit replaced by R, one sees that the result will begin to significantly decrease only when the 
position of the first node of the last Bessel function J m - mi (|k2|r) falls outside of the integration 
range, that is for |k 2 | < \m — mi\/R, where m x ^ m. If m x = m, then at |k 2 | <C l/R the last 
Bessel function can be approximated by the unity. Therefore, the limit |k 2 | — > taken in ( I73|) 
implies that 

| k 2 1 — > and R — > oo provided that | kg | R «C 1 ■ 

If instead R — > oo at fixed |k 2 |, then the transition of non- forward to forward results is 
discontinuous at |k 2 | = 0. 

In [5] it is claimed with the specific example of the Compton cross-section that if the 
transverse momentum transfer is small compared to k (in our notation, finite |k 2 | <C k at 
infinite R), then the mi-dependence has a narrow distribution peaked at mi = m. Our 
analysis does not support this claim. 

Looking back at the formalism used, we can conclude that our result that all mi essentially 
contribute to the decay rate/cross section just reflects the unfortunate choice of the same 
common axis z for all the twisted states appearing in the process. It does not give a clue of 
how twisted the final particles are with respect to their own propagation axes defined by their 
average values of the 3-momentum operator. Indeed, even a simple non-forward plane wave 
when expanded in the basis of twisted states contains all partial waves, see ([7]). Nevertheless 
it carries a zero orbital angular momentum with respect to its own direction of propagation. 
Therefore, it appears that a more physically reasonable quantity is the "orbital helicity", 
projection of orbital angular momentum on the axis of motion. The relevant question is then 
how this "orbital helicity" , not the OAM with respect to a fixed axis, is transferred from the 
initial to the final twisted state. We postpone this question for future studies. 

6.3 Conclusions 

Orbital angular momentum (OAM) is a new degree of freedom, which can be used in high- 
energy physics to gain more insight into properties of particles and their interactions. In this 
paper, focusing on the scalar case, we studied high-energy collisions in which initial and/or 
final particles were described by twisted states, i.e. carried a non-zero OAM. We derived 
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expressions for the cross sections for the single-twisted and double-twisted cases as well as for 
the two-particle decay rate of a twisted scalar, and observed a number of remarkable features. 

The single-twisted differential cross section is represented via the plane wave cross section 
averaged over the azimuthal angle of one of the incoming particles. If the initial twisted 
particle is prepared in a superposition state with different orbital quantum numbers, a Fourier 
analysis of the plane wave cross section can be performed. The expression we found for the 
double-twisted cross section is more intriguing, as it involves not only the plane wave cross 
section, but also the autocorrelation function of the amplitude. We stress that these features 
do not rely on large conical momentum spread in the twisted states, and their experimental 
study looks feasible even with today's technology. 

When analyzing the decay rate of a twisted particle, we established the procedure that 
allows one to pass from the plane wave to the twisted particle basis for the final state. These 
results can be now used, for example, to investigate the Compton backscattering of the twisted 
photons in the non- forward region, which was missing in the original suggestion [3 |6]. 

Analyzing parameters of the final twisted state, we discovered that there is no smooth 
transition from the non-forward to forward scattering. We attribute this result to the infinite 
transverse size of twisted states and to the deficiency of our formalism in which all the twisted 
states are defined with respect to the same z-axis. We expect that the problem will disappear 
if the orbital angular momentum projection is calculated not on the fixed reaction axis but on 
the direction of the outgoing twisted particle. Incorporation of this "orbital helicity" into the 
present formalism yet remains to be done. 
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A Restricted angular integrals 

Here we derive some properties of the restricted angular integrals which appear in the paper. 
The basic integral has the following form 

J = d<M<M (2) (ki + k 2 -k) (80) 
Jo 

and corresponds to a situation when two transverse momenta k x and k 2 of fixed absolute 
values sum up to a fixed momentum k with modulus |k| = k and azimuthal angle 

Clearly the integral can be non-zero only if it is possible at all to form a triangle with sides 
k, |k|i, and |k 2 |, that is, when |ki| and |k 2 | satisfy the "triangle rules": 

k < |ki| + |k 2 | , |k x | < «+ |k 2 | , |k 2 |<«+|ki|. (81) 

This allowed region on the (|ki|, |k 2 |)-plane has the form of a stripe shown in Fig. [2j If all 
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Figure 2: The allowed kinematical region of the values of |kx| and |k 2 | for a fixed k defined 
by the "triangle rules" (1ST]) . 

moduli are fixed, the integral flHUl) receives contributions only from two points (see Fig. [T] in 
the main text). Let us introduce 

S 1 = arccos ( - k2 ) , 5 2 = arccos ( - ^ ^ kl ) , (82) 

\ 2k|ki| / \ 2/t|k 2 | / 

Then, the two configurations of transverse momenta correspond to 

<\>k - 01 = ±<5i , 0fc - 02 = =F^2 , (83) 

so that the signs of 0^ — 0i and 0& — 02 are always opposite. 
To evaluate the integral itself, we rewrite the delta-function as 

<j(a)( kl _ ( k _ k2 )) = 2 5[kj - (k - k 2 ) 2 ] 5(0! - k _ k2 ) . (84) 
The 0i integration is then eliminated, and we obtain 

Jo = 2 J dcj) 2 5 [kj -K 2 -k 2 + 2n\k 2 \ cos(0 2 - fc )] = ^ , (85) 

where 



A = - yj 2(kf k 2 2 + K 2 k 2 + K 2 k2) _ K 4 _ k 4 _ k 4 (gg) 

is the area of the triangle with sides k, \ki\ and |k 2 |. When needed, this area can be also 
rewritten as 

A = i|ki||k 2 |sin(5i + 5 2 ). (87) 

Let us now consider an integral similar to (|80|) but with an extra angular dependence in 
the integrand: 

"2n 



j= r d<h # 2 6^ ^ + k 2 - k) • /(0!, , 

JO 



21 



with a regular function /(0 1; 2 ). Here again the momenta must satisfy the triangle rules 
(jSI) and the integral receives contribution only from two points (183]) in the (<pi, 2 )-space. 
Denoting the values of f(4>i, 02 ) & t these two points as /+ and /_, respectively, we obtain 

J=^. (89) 



We can also rewrite this result as 

c2w 



J= f #i# 2 £ (2) (ki + k 2 -k) ■ 



(90) 



which means that in these circumstances the delta-function plays the role of a functional that 
maps a test function f((f>x, 2 ) to its symmetrized value (/ + + /_)/2. 
In the particular case of / = exp(im 1 1 + 2m 2 2 ), (18"9"j) take form 

j = e i(m 1+m2 )4> k cosjmA - m 2 £ 2 ) ^ (gi) 

B Regular izat ion of \I\ 2 

Here we calculate the large-i? behavior of the mi-sum of the squares of the triple-Bessel 
integral: 



(92) 



which appears in the decay rate (|74l) . Evaluation of the integral itself with R — > oo performed 
in the main text shows that it can be non-zero only if k, k±, K2 satisfy the triangle rules (181 p . 
i.e. a triangle with these sides can be constructed. Since m describes the initial state, we take 
it small and not growing with R: m <C m\ max = n±R, while mi can extend up to vn\ max . The 
final Kx, k>2 can be much larger than k. 

Since the expression (1921) is regularized with large but finite R, the summation and inte- 
gration can be interchanged: 

rdrr'dr' J m (Kr)J m (Kr') ^ J mi (Kir) J m - mi (K 2 r) J mi (Kir') J m - mi («2r') . (93) 

Thanks to the properties of the Bessel functions, only mi's up to mm(Ki j2 r, K\^r') are ef- 
fectively contributing to this sum; for larger m\ the Bessel functions strongly decrease. But 
r, r' < R, which means that the limits on the summation can in fact be safely extended to the 
infinity. Then, the sum of the product of four Bessel functions is treated in the following way: 

+ 0O 

^2 4j )4-mi («2'")4 1 4-mi («20 (94) 

+ 0O 

u mi ,m. 



mi =— oo 

+oo 



mi ,m,=— oo 
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ZJt_^ 271 Jo 

+00 



mi ,m, =— 00 



£ e fmia J mi (Kir)J m _ mi (K 2 r) 

mi=- 00 



mi =—00 



The first sum in the square brackets is calculates as follows: 

+00 

e imia J mi (K ir )J m _ mi (K 2 r) 

m\=— 00 



(2tt) 2 

(-i) m 
2vr 



+00 

^ e iKircos</>i+ire 2 rcos<fe ^ ^ gimia+imi^i+ifm-mi)^ 

mx=— 00 



^ gjm0i gj«ir cos <f>i+ik2'r cos((/>i+a) 



(95) 



The combination of angles and momenta inside the exponential can be expressed as 

Ki cos 0i + k 2 cos(0i + a) = I ki + k 2 | Q cos(0i + 50) , (96) 

where 



ki + k2| a = y k>i + «2 + 2/«i/«2COS a , tan< 



^2 sin a 



(97) 



Ki + k 2 cos a 

Geometrically, |ki + k2| Q is the norm of the sum of two vectors of moduli K\ and k 2 and the 
relative azimuthal angle a. Therefore, (|95l) is 



(98) 



e m ^J mi (K ir )J m _ mi (K 2 r) = e *»<«-W J m (|ki + k 2 | a r) . 

rrti=— oo 



This expression can be viewed as a 2D generalization of the well-known addition formula for 
the Bessel functions 

+00 

mi=— 00 

Now, the second sum differs only by a — > —a (or, alternatively, complex conjugation) and 
r — > r'. Therefore, the summation (1941) is simplified to 



1 f 27T 

— / da J m (|ki + k 2 | a r) J m (|ki + k 2 | a r') . 
^ Jo 



(99) 



Note that this integral involves only the Bessel functions of small order. We now plug this 
expression in fl93l) and get 



-2n 



27T 



da 



rdr J m (nr) J m (|k x + k 2 | Q r) 



J r'dr' J TO («r') J m (|ki + k 2 | a r') 



(100) 
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As usual, we extend the integration range in one of the integrals to infinity, which gives a 
delta-function, and then we use it on the second integral calculated up to R. The original 
expression fl92l) then becomes 
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(101) 



where A is, as always, the area of the triangle with sides n, k±, k^. 
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